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Abstract. We study the quantum forces that act between two nearby 
conductors due to electronic tunneling. We derive an expression for these forces by 
calculating the flux of momentum arising from the overlap of evanescent electronic 
fields. Our result is written in terms of the electronic reflection amplitudes 
of the conductors and it has the same structure as Lifshitz's formula for the 
electromagnetically mediated Casimir forces. We evaluate the tunneling force 
between two semiinfinite conductors and between two thin films separated by an 
insulating gap. We discuss some applications of our results. 
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1. Introduction 

The increased accuracy of experimental studies 012 El EI H] of the Casimir force [H] 
between conducting bodies has opened the possibility of exploring new ideas related 
to the understanding and control of quantum vacuum fluctuations. Research projects 
on the dynamical Casimir effect ,7. , Casimir torques [5] , or the possible applications of 
the Casimir forces in the development of micro and nano-electromechanical devices are 
now under way Understandig Casimir forces has become fundamental in the 

investigation of deviations of Newton's Gravitational Law at micrometer dimensions, 
related to the search for extra dimensions in space-time . 

The usual Casimir effect may be understood as a force due to the quantum nature 
of the electromagnetic radiation. In this paper we study another source of quantum 
forces, namely, the tunnel effect. Particles that are able to tunnel across a barrier have 
more space available to them. Thus, their contribution to the total energy of a many 
body system such as two metallic slabs with neighboring surfaces may diminish. As 
the amount of space gained and the number of particles capable of tunneling depends 
on the width of the barrier, besides its height, there must be a force that performs work 
when the width is modified. Since quantum tunneling arises from evanescent electronic 
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fields, this force is similar to the contributions of evanescent electromagnetic waves to 
the standard Casimir force. 

In this paper we concentrate our attention on conduction electrons in conductors, 
that is, on massive non-relativistic Fermions. We derive the tunneling force by 
calculating the flux of momentum between two regions delimited by an arbitrary 
potential V(x). We express the momentum flux in terms of the Green's function of 
the system, which we evaluate by means of a scattering method involving amplitude 
reflection coefficients ^JEI|- This method yields an expression for the tunneling force 
with a structure that is essentially identical to Lifshitz's formula ^3]- We first perform 
the calculation for a one-dimensional system. We then extend the calculation to the 
three-dimensional case. Finally, we evaluate the tunneling force for a configuration 
consisting of two semiinfinite or two thin metallic slabs separated by a thin insulating 
gap and we discuss some applications of our results. 



2. One dimensional systems 

The dynamical equation for the wavefunctions of an electronic system may be derived 
from a Lagrangian density, |lf>| 

c = + vm 2 + y (we* - r^t)) , (1) 

for which Euler-Lagrange's equations yields Schrodinger's equation, 
„ dC dC dC _ , h 2 , 

dt Wn) + dz m~) ' w = lhtt + - v * = °- (2) 

The wavefunction carries mechanical properties which may be derived from C: we 
may obtain a momentum density 

o-y-'W^W*- -?<**! "*•*->■ < 3 > 

where c is the speed of light in vacuum, as well as a momentum flux |16j 

t ' = + - c - £ fc|2 - m * - " ( ***' - ♦•**) ■ (4) 

Consider now an eigenstate ip n of the Schrodinger Hamiltonian H = 
— {h 2 /2m)&l + V corresponding to an energy E n . From Eq. its contribution 

to the momentum flux is 

T z * = ^\d z Tp n \ 2 + (E n -V)\i> n \ 2 . (5) 
2m 

Within a region V in which V{z) may be taken as a constant, we can write E n — V = 
h 2 kf l /2m, i.e., the kinetic energy of particles with wavenumber k within V. We now 
sum the contributions (JSJ over all the occupied orbitals, 

T ^ Z ) = ^E / dE6(E-E n )f(E)[k 2 \Mz)\ 2 + \dMz)\ 2 }, (6) 
Ira L — ' / 
n J 

where f(E n ) is the occupation number of orbital n, given in equilibrium by the Fcrmi- 
Dirac distribution function, and k 2 = 2m(E — V)/h 2 . The energy integration and 
Dirac's 5 allow us to write T* in terms of the Green's function of the system, 



G E (z, z') = (z 



E-H' 
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employing the relation Im(E + — E n ) 1 = —7r5(E — E n ), where E + = E + ir\ with E 
and i] — > + real. Substituting this latter relation in © and employing (JJJ we obtain 

T:(z) = -^Im / d£;[fc 2 G is+ (z,z') + ^9 Z 'G £+ (z,z')] 2 ^ 2 /(^)- (8) 
Notice that Eq.© may be interpreted as 

T*{z) = J dE P % f f(E)(hk) (^j , (9) 

where ±hk is the momentum of a particle which moves at velocity ±hk/m, thus 
contributing the amount (hk)hk / m to the momentum flux, and 



Pe( z ) = -^~ lm 

Z7T 



G E (z,z') + -^d z d z >G E {z,z') 



(10) 



plays the role of an effective local density of states. 

We now assume that V has a width L and is bounded on both sides by arbitrary 
potentials, and we evaluate the Green's function following a scattering approach 
PI ITS] , Within V, the solution of (E - H)G E (z,z') = 5{z - z') may be written 
as Ge(z,z') = (2m/h 2 )'ijjL(zL)^pR(zii)/W, where ipL and ipn are the two solutions of 
the Schrodinger-like homogeneous equation (E — H)ip = that satisfy the boundary 
conditions on the left and the right side of the system respectively, W = ipLipR.z — 
ipL,zipR is their Wronskian, and z E and Zr are the smallest and the largest among z 
and z'. We write ip L (z) = e~ ikzh + r x e ikZL and ip R (z) = e ik( - ZR - L ^ + r 2 e- ik ( ZR - L \ 
where r\ and r 2 are the reflection amplitudes for particles impinging on the left and 
right boundaries of V, which we assume at z = and z = L, and we obtain 

, 2m (e~ ikZL + r ie tkZL ) ( e lk( - ZR - L + r 2 e- ik ( ZR - L ) 

Ge(z,z) = —2 — _ ikL , 2ikT\ ' (•'■•'•) 

% 2ike %Klj (1 — r\r 2 e ) 

which together with Eq. (JSJl yields the momentum flowing within V, 

iRc [ d Ek 1 + r ^ e T L 
7r J 1 - rir 2 e z 

which may be conveniently written as 

■Km J 1 — r\r 2 e l 

by using the relationship E — Ti 2 k 2 /2m + V and changing integration variable. Notice 
that, as expected in an equilibrium situation, T z is independent of z. 



Tt = -Re / dEk J^^ 2ikL f(E), (12) 



T Z Z = / dkk 2 J ± \ 2lkL f(E) (13) 



3. Three dimensional systems 

The generalization of the results derived above to the three dimensional case is 
straightforward for systems which are translationally invariant along a symmetry 
plane, say xy. In that case, the parallel wave vector Q — (Q x ,Qy) is a conserved 
quantity, and for each Q the problem is identical to the ID case. Thus, we only have 
to sum Eq. \12\ over the allowed wavevectors, 

V = T^- Re / d 2 Q { dkk 2 ] + rir2e l lk k L r f(E) (14) 
4n 3 m J J \-nr 2 e 2lkL 
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where we introduced the number Ad 2 Q/(2ir) 2 of wavevectors within a region d 2 Q of 
reciprocal space by applying Born- von Karman boundary conditions in a system with 
total area A — ► oo, and we introduced the momentum flux density T z z = T* / A. Notice 
that —T" coincides with the zz component of the stress tensor as defined in elasticity 
theory. 

We remark that the structure of Eq. 114JI is essentially identical to Lifshitz's 
formula for the Casimir effect between two materials when written in terms of 
their optical reflection coefficients EJ. The main differences are that the 
electromagnetic field has two independent transverse polarizations whose contributions 
would have to be summed over, and that the speed of light is a constant c, while 
the speed of electrons is proportional to the wavevector, i.e., the dispersion relation 
between electrons and photons are different, and consequently, there is an extra power 
offrinEq. ifTI]!. 

As r\ and r-2 are independent of Q for scalar fields, the first integral in Eq. (|14|l 
may be performed immediately. At zero temperature we obtain 

V = -I Re f dk U F - fk 2 ) e ) + rir f^ L , (15) 
7r z J \ 2m J 1 — rir2e ZIKIj 

where the integration region includes all states below the Fermi level, whose kinetic 
energy within V is Kp, and for which we took f(E) — 2, including the spin degeneracy. 



4. Applications 

4--1- One semiinfinite metal 

Within the bulk of a semiinfinite metal the electrons are reflected by the surface 
potential barrier on one side, while there is no barrier on the other side. Thus, the 
pressure p with which the electrons push the surface of the metal may be obtained by 
setting r% = in Eq. (|15fl . The result is simply 

1 r kp f h 2 k 2 \ 9 

P = V= ^l dk [ Ef ^ = 5 F > (16) 

where hkp is the Fermi momentum, Ep — Kp (within the metal) is the Fermi energy, 
and n = kp/Zir 2 is the electronic density. As could have been expected, this result 
coincides with the well known pressure of a degenerate fermion gas |17) . 



4-2. Two semiinfinite metals 

We consider now two identical semiinfinite metals separated by vacuum. The force 
F/A per unit area between both metals may be obtained from the momentum flux 
(1151) within the vacuum region, where the wavefunction of all the occupied states are 
evanescent, and it may be written as 

^ = -2Im^/J^ ( 4-.V^p (17) 

where = r 2 e ~ 2KL , and we wrote the wavenumber k = in in terms of the decay 
constant k. The integration limits in (|17|l are the decay constants for electrons at 
the bottom of the conduction band, kq = y^[2m(W + Ep)/h 2 \, and at the Fermi 
level kf = \f(2mW/h 2 ), while W = Kp (within vacuum) is the work function, 
and r — r\ — r2 is the complex reflection amplitude corresponding to evanescent 



(M 
-Si 



Casimir-like tunneling-induced electronic forces 
10 
1 

0.1 



0.01 
0.001 
le-04 
le-05 b- 



le-06 



; 1 1 1 1 


1 1 1 1 


i : 


\w = w '""---.^ 


ocF x (10) 


- 


~-w = o.i ^v;;^ 






=" "cx Fi(O.l) " 






I I I I 


I I I I 


1 



0.1 0.2 0.3 0.4 0.5 0.6 

Distance kpL 



0.7 0.8 0.9 



Figure 1. Force por unit area between two scmiinfinite metals as function of their 
separation L for different values 0.1, 0.3, 1, 3, 10 of the dimensionless work function 
W = W/2Ep. Also shown schematically is the contribution e~ 2KpL <x F\ 
expected for one electron at the Fermi surface in the cases W = 0.1, 10. 



wavefunctions that propagate (i.e., decay) through vacuum towards a surface and are 
reflected back. Assuming that the potential V(z) is constant within the metals and 
within vacuum, and that it changes abruptly at the vacuum-metal interface by an 
amount W+Ep, the reflection amplitude may be calculated as r — (in— A;jvf)/(i/t+fcM), 
where ku = yJ[2m(W + E p) / ft 2 — n 2 ] is the wavenumber within the metal of the state 
corresponding to k. 

In Fig. ^we plot the force per unit area as a function of distance for different 
values of the workfunction W. The force is attractive, seems to decay exponentially 
for large separations and attains a finite value at zero separation. For large W the 
force is larger at small separation and smaller at large separations as the energy 
decays very fast towards that of two isolated semiinfinite metals. As could have been 
expected, the smaller the work function, the larger the spatial range of the force. 
We might expect the decay to be dominated by those electrons closest to the Fermi 
energy whose contribution becomes proportional to e ~ 2tiFL . Fig. ^includes two curves 
illustrating this behavior for the cases of large and small W. The actual decay of the 
force is slightly faster, more so for small W . This is due to the fact that not only 
the contribution of each electron decays with increasing distance, but also the number 
of electrons that contribute effectively to the force. Furthermore, the phase space 
available right at the Fermi energy is null, due to the prefactor k 2 f — n 2 in Eq. i|17|) . 
so the contributing electrons have a slightly larger decay constant (i.e., smaller range) 
than those at the Fermi level. 

In Fig. we show the force for several distances as a function of W . For finite 
separation distances the force is small when W is large, as the surfaces don't feel each 
other anymore, and it is also small for small workfunction, as the electrons tunnel too 
easily between the two metals, not caring about the separation. Thus, the magnitude 
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Figure 2. Force between two metals as a function of the work function for several 
values of the separation L = 0, 0.1, 0.2 in units of k~ . 



(18) 



of the force is largest at some intermediate value of W which increases as L decreases. 
At contact, L = 0, there is no such extreme anymore and we obtain a linear behavior, 

F(0) = _g| (Ep W" 
A ir 2 \ 5 : 

as can be shown by integrating Eq. (|17f) analytically. F(0) is the force that would 
be required in order to break an infinite metal into two semiinfinite ones. Eq. (|18fl 
actually overestimates the ultimate breaking strength of real materials by several 
orders of magnitude as our model fails to account for dislocations whose motion 
within the metal would relax the stress, and for the growth of fractures which are 
actually responsible for the failure of real metals. Real metals break gradually, not 
simultaneously over the whole separation surface. Nevertheless, integrating Eq. (|17|l 
over L we have obtained an analytical estimate of the surface energy of metals in terms 
only of their Fermi energy and their work function. This turns out to be surprisingly 
accurate |18] given our simplifying assumptions, namely, our use of an independent 
free particle model, neglecting the crystalline structure, the electronic charge, and 
many body corrections, as well as our use of a square potential barrier at the surface. 



5. Thin films 



Eq. (|17|) may be employed to calculate the force between more complicated systems 
simply by introducing the appropriate value of the reflection amplitude. For example, 
in Fig. 0|we display the force between two free standing very thin metallic films as a 
function of distance for a fixed Fermi energy, or more properly, a fixed electrochemical 
potential. Notice that for very thin films the force is identically zero, as there are no 
states at all below the Fermi energy and therefore there are no available electrons to 
tunnel between the films. For wider films the force is finite at small separations but 
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Figure 3. Force between two films of widths d = 0.5, 1, 2, 4, 
as a function of the distance L between them. We took W = 
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becomes zero after a finite separation. As the width is further increased, the force 
approaches that corresponding to semiinfinite metals, although not monotonously; it 
actually oscillates between larger and smaller values. This behavior may be understood 
by considering the finite size effect on the levels of the individual films, and the 
interaction of levels within both films, yielding alternating bonding and antibonding 
states which may be occupied only when they lie below the Fermi level. 

6. Conclusions 

By calculating the mechanical properties carried by the electronic wavefunctions, we 
have shown that the interchange of electrons between conductors produces a force 
that may be calculated in terms of the electronic reflection amplitudes using formulae 
that are very closely related to Lifshitz formula for the usual Casimir force. We 
illustrated our formalism by calculating the electronic pressure within a conductor 
and the force between semiinfinite conductors and between thin films at very small 
distances, of the order of the Fermi wavelength. These distances are extremely small, 
beyond the expected limit of validity of the usual Casimir effect. Thus, we expect 
our results to be important to study the forces that act, for example, between the 
tip and the substrate of a scanning tunneling microscope |19j . We discussed how 
our results may be employed to calculate the surface energy of conductors without 
having to substract total energies. Other applications which are currently under study 
include the calculation of the force between impurities embedded within three and one- 
dimensional conductors. Although our current calculations were performed for free, 
independent electron conductors, we believe that our scattering approach might be 
generalized to more realistic systems of interacting electrons. 
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